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Abstract 

The gluon polarization AG/G is the key to a further clarification of the spin 
structure of the nucleon. The COMPASS collaboration at CERN has set out to 
CN ■ undertake the direct determination of this quantity. It accesses the gluon distri- 

ct . bution via the photon-gluon fusion process (PGF) in scattering polarized muons 

off a polarized deuteron target. And it explores three different channels to tag the 
PGF: open charm production and high transverse momentum (high pr), in either 
electroproduction (Q 2 > 1 GeV 2 ) or quasi-real photoproduction (Q 2 < 1 GeV 2 ). 

The high pj- quasi-real photoproduction channel yields the most precise mea- 
surement. The result indicates that the gluon polarization is small. I describe 
its experimental aspects and its theoretical framework, based on PYTHIA. And I 
report on the preliminary results obtained in the other two channels, and on the 
prospects for future analysis and data taking. 

o ' 
o 



1 Introduction 



<D The spin 1/2 of the nucleon can be decomposed as follows 

- = -AZ + AG + L Z 
X : 2 2 

where the right hand side terms designate the contributions of the spin of the quarks, the 
spin of the gluons and the angular momentum of the quarks and gluons, respectively. 

In the recent years much effort has been put in determining AS. This quantity can 
be derived from the measurement of the spin dependent structure function gi by polar- 
ized inclusive deep inelastic lepton-nucleon scattering {DIS) experiments. Measurements 
were carried out at CERN, SLAC, DESY and JLab. COMPASS itself is presently car- 
rying on with this program. The results lead to the conclusion that AE is surprisingly 
small, significantly smaller than predicted by the Ellis- Jaffe sum rule [1] for example. An 
evaluation of AS taking into account the latest data points measured by COMPASS is 
presented by R. Windmolders at the SPIN05 workshop [2]. 

A solution to the problem was put forward in 1988 [3, 4, 5]. It involves a leading 
order contribution to the polarized DIS cross-section originating from the axial anomaly 
of QCD, a s /2ir AG, which is anomalous in the sense that it does not vanish in the 
asymptotic limit: in leading order evolution AG grows with InQ 2 whereas a s is inversely 
proportional to InQ 2 . This anomalous gluon contribution introduces some freedom in 
the definition of AS, but it can in any case reconcile polarized DIS data with QCD 
predictions given a large enough, positive, value for the first moment of AG [6]. 
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The unpolarized gluon distribution, G, can be determined from the dependence of the 
inclusive DIS cross-section upon Q 2 . In the polarized case, however, DIS data cover too 
small a range in Q 2 for this method to significantly constrain AG. 

A direct measurement of the gluon polariza- 
tion, AG/G, is therefore the most promising 
way to further clarify the nucleon's spin puz- 
zle. Several experimental projects have recently 
been started to carry it out, using different ap- 
proaches. In COMPASS, we access the gluon 
distribution via the photon-gluon fusion (PGF) 
process (depicted on Fig. 1), whereby a virtual 
photon couples to a gluon via a qq pair. And 
we consider two different selections of the PGF: 
open charm production and high hadron pro- 
duction. The two share a set of common fea- 
tures. Both have been successfully used to di- 
rectly measure the unpolarized gluon distribu- 
tion at the HERA collider experiments (with the 
difference that high p? refers to the production of 
jets there). And for both, the QCD scale can be 
set irrespective of Q 2 , by the charm mass and the 
Pt cut, respectively. But the two selections rep- 
resent diametrically opposed trade-offs between 
the conflicting requirements of statistics and pu- 
rity. Open charm is the purest. It provides a model-independent access to AG/G and 
for this reason remains our golden channel. I will present it first. And I will present next 
the high px case, which we subdivide into two sub-cases depending upon the Q 2 of the 
exchanged photon. I start with some experimental essentials. 




Figure 1. Photon-gluon fusion process 
in a muon nucleon scattering event. 



2 Experimental essentials 

The COMPASS spectrometer is described in details in [7]. 

Its experimental setup was designed to allow a precise determination of asymmetries. 
An important point in this respect, is the control of fake asymmetries. We achieve it 
thanks to the simultaneous measurement of both orientations of target spin in two oppo- 
sitely polarized target cells, u and d, and to a frequent reversal of spin orientations, so 
that fluctuations in acceptance and incident muon flux cancel out in the formula for the 
counting asymmetry A: 

2 [ N r + ivf + ivtv ' [ ) 

where f|~f and "fiT denote the two spin configurations. (Note that weighted asymmetries 
are used instead of (1) in all calculations presented below.) 

The cross-section helicity asymmetry, A\\, is related to the counting asymmetry by 
factors describing the polarization of the incoming particles, Pb for the beam (~76%), Pt 
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and / for the target polarization (~50%) and for the, process dependent, dilution factor 
(~40%). It is best expressed as A\\/D, 

A\\/D = A/ (P b xP t fxD) , (2) 

where one takes also into account a kinematical factor, D, describing the polarization 
transfer from the muon to the photon. D is process dependent and typically averages to 
-60%. 

During its three years of running from 2002 to 2004, the experiment has accumulated 
~2.4 fb^ 1 of data with its target polarized longitudinally. The results presented here 
correspond to about half of these data. 



3 Open Charm 

This channel was discussed by many authors [8, 9] as a good candidate to access AG/G. 
Since there is no or only a small intrinsic charm in the nucleon, diagrams with an incoming 
charm quark do not contribute and PGF enters alone in leading order. 

In COMPASS we tag open charm, and hence PGF, by the reconstruction of a D° 
meson. The D° — > Kir decay channel is used, with K and it identified in the RICH. The 
main difficulty lies in the associated combinatorial background. This is a major concern 
in our experiment, where the vertex resolution is not sufficient to resolve the decay vertex 
from the primary vertex, because of the thickness of the target. 

Special care is therefore taken to optimize the use of the data. First, the favorable 
cases when the D° comes from a D* — > D°n decay are counted separately, cf. Fig. 2. 
Secondly, a weighting procedure is applied for the derivation of AG/G: 

AG/G ~ T T EfV + EfV 1 ~ (i + B/S h (3) 

where an is PGF's analyzing power and S and B are the signal and background counts. 

A Monte Carlo simulation of the experiment is used to calculate estimates of all 
quantities that cannot be directly computed from the hadron level kinematics, such as 
a LL and x g , the momentum fraction carried by the gluon. It is based on the event generator 
AROMA [10], which uses the PGF matrix element to generate charmed hadrons. an is 
first parameterized as a function of the transverse momentum and energy fraction of the 
reconstructed charmed meson, and of a kinematical factor, D, describing the polarization 
transfer to the virtual photon, so that it can be estimated on an event by event basis in 
equation (3). 

In these conditions, a preliminary analysis, bearing on ~l/2 of the data, yields: 

AG/G = -1.08 ± 0.73(staf .) at x g = 0.15 ± 0.08 RMS. 

This means a clear lack of precision, even when one takes into account all of the 
available statistics. Obviously more data are needed. COMPASS is going to resume data 
taking in 2006, with an improved experimental apparatus. Several upgrades, targeting 
three important aspects of the D° selection will be introduced: 
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D* candidates 2003 




Figure 2 D° peak in the Kir invariant mass distribution for all events (top) and .D*-tagged 
events (bottom) for 2003 (i.e. ~l/3 of total) data. S e ff is the effective number of events 
S e ff = S 2 1 (S + B) where S and B are signal and background counts. A bump shows up at low 
mass in the D* case, attributable to D° — > Kirir°. It is not included in the S count. 



• An enlarged angular acceptance, thanks to a new target superconducting solenoid 
(with an aperture of 180 mrad instead of 70 mrad), and a better coverage of large 
angles by drift chambers. 

• A better particle ID, thanks to a faster RICH detector, able to clear away the high 
level of background generated by the halo of our muon beam. 

• Electromagnetic calorimetry, that will allow us to explore D° decays with a tt° in 
the final state, mainly D° — > Ktttt° with a BR of 13% against 3.8% for D° — > Kir. 



4 High pt 

The alternative approach to select PGF is to require hadron production with a high 
transverse momentum with respect to the virtual photon [11]. This suppresses leading 
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order 7*g — > q events, where the fragmenting quark goes into the direction of the photon. 
However the suppression is not perfect. And other competing second order processes have 
to be taken into account. These include the QCD Compton process, and, mainly at low 
Q 2 , processes involving resolved photons. In these conditions, and limiting oneself to 
leading order, where the different processes are well separated, the cross-section helicity 
asymmetry reads: 

A\\ = RpQ F a LL AG/G + A 

background j (4) 

where Rpgf is the fraction of PGF events, a I [ ( [ F is PGF analyzing power and Ab ac kground 
describes the contribution of the background processes. In order to retrieve AG/G, Rpgf, 
Abackground and the parton level kinematics defining a P( [ F need be determined by a simula- 
tion of the experiment. In COMPASS, we resort to the event generators LEPTO [12] and 
PYTHIA [13] to do the computation, which introduces a model dependence in our mea- 
surement. LEPTO covers the electroproduction regime at Q 2 > 1 GeV 2 and PYTHIA 
the quasi- real photoproduction regime at Q 2 < 1 GeV 2 . In both cases, the leading order 
approximation is consistently applied, in the choice of PDFs (we use GRV98 [14] and 
GRSV2000 [15] at leading order), in the calculation of the analyzing powers and by turn- 
ing off parton showers in JETSET. And for both, the event selection follows a same path. 
Two hadrons are required, with px > 0.7 GeV and Epf- > 2.5 GeV 2 . Contributions 
from target fragmentation and resonances are excluded by requiring xf > 0.1 and the 
invariant mass of the two-hadron system to be larger than 1.5 GeV. But the two cases 
are attractive in their own right: Q 2 < 1 GeV 2 yields much higher statistics, a factor 10, 
but it also involves more model dependence. 

4.1 High p T , Q 2 < 1 GeV 2 

The details of the measurement are available in [16]. I recall here the main aspects. 

PYTHIA provides a complete coverage of photon-nucleon interactions at low Q 2 [18]. 
We checked that it reproduces reasonably well our own data, by propagating generated 
events through a GEANT model of the spectrometer and processing them through the 
full reconstruction and analysis chain. A good agreement was readily obtained for the 
kinematical variables Q 2 and y. A single parameter, the width of the intrinsic kj* of 
the partons in the photon, had to be adjusted (from 1 to 0.5 GeV) to fit the longitudi- 
nal and transverse momentum distributions of the two hadrons. As an example, Fig. 3 
shows the agreement obtained for the pr distribution of one of the hadrons. The relative 
contributions, R, of the various processes thus generated are shown on Fig. 4. 

To compute the asymmetry of a given process, it is mandatory to find a hard scale 
allowing the factorization of the cross-section into a hard partonic cross-section calculable 
perturbatively and a soft parton distribution function. In our case, this hard scale is to 
be looked for in the pr of the partons. Thanks to the hadron p T cut, for most events, 
one parton at least turns out to fulfill the condition. However this is not true for leading 
order 7*g — > q events, and for some other events which PYTHIA classifies as "Low-p-r" • 
In these conditions, we can expand equation (4) as follows: 



5 




Figure 3 Comparison between data and Monte Carlo. The pt distributions for the hadron with 
highest pt are displayed, superimposed, for the two of our triggers. The bottom plots show the 
ratios: data divided by Monte Carlo. 



f,f~f=u,d,s,u,d,s,G 



- tiPGF { D ) -g- + tiQCDC { Y) Al ' + 

^// 7 i^LL j j^y ) -^Leading AL eac n n g ~\~ Rlow—pt Al w—pt > (p) 



where the QCDC term makes use of the inclusive photon-deuteron asymmetry, A\, and 
the resolved photon term makes use of the PDFs in the deuteron, / and A/, and in 
the photon, f 1 and A/ 7 , and where D is a kinematical factor which approximates the 
amount of polarization transfered from the muon to the photon. The last two terms, 
corresponding to the leading and low pr processes, cannot be expanded, since there is no 
hard scale allowing to factorize their cross-section factor. However, the asymmetry for 
this kind of events is known to be small, as indicated by measurements at low Q 2 [17], and 
they account for only 7% of the high pt sample, cf. Fig. 4. Therefore, they are neglected. 
All the PDFs appearing in the remaining terms, apart from the targeted AG, have been 
measured, except the polarized distributions of the photon A/ 7 . However, each of these 
is bound by a minimal and a maximal saturation scenarios, as described in [19] 1 . We can 
therefore solve equation (5) to determine AG/G from the measurement of (An/D). Note 
that the PGF term is not the only one containing AG/G, the resolved photon term is 

1 Our definition of the minimal scenario corresponds to a maximally negative polarization of the parton 
in the VMD part of the photon. 
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Figure 4 Relative contributions R of the dominant PYTHIA processes to the Monte Carlo 
sample of high pj* events at Q 2 < 1 GeV 2 . Left: direct processes, right: resolved photon 
processes. Longitudinal photons, as well as some minor resolved photon contributors, are not 
shown. 

also doing so. The two contributions have opposite signs and cancel each other to some 
extent: the PGF term is equal to —0.29 x AG/G and the linear coefficient of AG/G 
in the resolved photon term is positive, ranging from +.012 to +0.078 depending upon 
which saturation scenarios are used. 

The Monte Carlo analysis outlined above introduces systematic uncertainties. We have 
considered the following two sources: NLO effects and the tuning of PYTHIA. We have 
estimated their order of magnitude by repeating the analysis several times with modified 
Monte Carlo parameters. For the NLO effects, the renormalization and factorization scales 
were multiplied and divided by two and the parton shower mechanism was activated. For 
the other source, we have explored the parameter space, while maintaining a reasonable 
agreement between simulated and real data. The most critical parameter appeared to be 
the width of the intrinsic kq* in the photon. 

An/D can be obtained from equations (1),(2). In order to minimize the statistical 
error, we use instead the weighted method described in [20] . With the high px selection 
defined above, and the added condition 0.35 < y < 0.9, we obtain: 

(A\\/D) = 0.002 ± 0M9(stat.) ± 0.003(syst.) , 

where the systematic uncertainty accounts for false asymmetries, which were estimated 
using a sample of low px events with much larger statistics. 

For the derivation of AG/G, we first consider the minimum and maximum values 
obtained from equation (5) with the various saturation scenarios for the polarized PDFs 
of the photon mentioned above: 



7 



( AG/G ) min = 0.016 ± 0.068(stat.) ± 0Ml(exp.syst.) ± 0M8(MCsyst.) 
( AG/G ) max = 0.031 ± 0.089(stat.) ± 0.0U(exp.syst.) ± 0.052(MCsyst.) 

This leads to the central value: 

AG/G = 0.024 ± 0.089(stat.) ± 0.057(syst.) at x g = 0.095^ and /j, 2 ~ 3 GeV 2 . 

where the difference between the min and max values is included in the systematics and 
all systematics are added quadratically. x g and the average scale fi 2 are both obtained 
from the Monte Carlo simulation. 

4.2 High p T , Q 2 > 1 GeV 2 

This time the analysis is only preliminary. It follows almost the same path as the one just 
described for the low Q 2 case. And I will only mention what singles it out. 

The Q 2 cut introduces several simplifications to equation (5). It provides a hard 
scale for all processes. This eliminates all purely soft processes previously classified as 
"Low-p r ". And it allows us to factorize the leading process term into a hard cross- 
section and a soft part, which we approximate by A\ as for the QCDC term. Besides it 
reduces drastically the cross-section of the resolved photon processes, which are altogether 
neglected. 

Further simplification is obtained by restricting the analysis to the low Bjorken x 
region, x B < 0.05, where A x for the deuteron is small [21]. Hence, the contributions of 
the QCDC and leading processes to the asymmetry become negligible. 

After all selection cuts, we obtain: 

(A\\/D) = -0.015 ± 0.080(stat.) ± 0M3(syst.) 

The Monte Carlo simulation is performed using LEPTO, with a set of fragmentations 
parameters modified to fit the data. It yields Rpgf — 34%, (a LL /D) ~ —75% and: 

AG/G = 0.06 ± 0M(stat.) ± 0.06(syst.) at x g = 0.15 ± 0.08 RMS , 

where the systematical error accounts only for the experimental uncertainty. 

5 Conclusion and outlook 

Our high p? results for the gluon polarization are plotted on Fig. 5. They are compared 
with previous direct measurements by SMC [22] and HERMES [23]. The SMC measure- 
ment uses high p T events at Q 2 > 1 GeV 2 , while the HERMES result is derived from 
data mostly at low Q 2 . Both neglect resolved photon processes. 

It is interesting to compare our most precise result, from high p? at low Q 2 , with 
distributions of AG/G(x g ) obtained from QCD fits. Three such distributions, from 
GRSV2000 [15], are plotted on Fig. 5. They correspond to three hypotheses on the gluon 
polarization at an input scale of \j? = 0.40 GeV 2 : maximal polarization (max), best fit 
(std) and zero polarization (0). The distributions are then evolved to fi 2 = 3 GeV 2 , 
where their first moments are 2.5, 0.6 and 0.2, respectively. Our low Q 2 result clearly 
favors parameterizations with a low gluon polarization. Two other recent NLO fits [24, 25] 
(not shown), with first moment values at /i 2 = 3 GeV 2 ranging from 0.26 to 0.8, were 
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Figure 5 Comparison of the AG/G measurements from COMPASS (the two high px results are 
shown, that at Q 2 > 1 GeV 2 being preliminary), SMC [22] and HERMES [23]. For each point, 
the vertical error bar includes the statistical uncertainty only and the horizontal bar represents 
the x g range. The curves show three parameterizations of AG /G{x g ) at [i 2 = 3 GeV 2 from the 
NLO fits of [15]. (see text for details). 



also checked against our result at x g = 0.095 and found to be compatible within 1.5 a. 
The anomaly mechanism mentioned in the introduction would require a first moment of 
AG of about 3 to explain the polarized DIS data. The small value of AG/G from our 
measurement cannot rule out this possibility. However it makes it very unlikely. 

The two high results discussed above assume the validity of PYTHIA and LEPTO, 
respectively. The open charm channel would allow instead to make a model independent 
measurement. It still suffers from a lack of statistics but remains our golden channel. 
COMPASS resumes data taking in 2006, with an improved experimental apparatus, and 
will strive to complete the measurement. 

The high pt channels will also benefit from the new data. They will have to include 
first the already existing data which remain to be analyzed. Which will bring down the 
statistical precision to 5 AG/G = 0.22 for the high Q 2 case and 5 AG/G = 0.06 for the 
low Q 2 case. For the latter, the aim will be to use the extra statistics brought by 2006 
data to create two bins in x g . 
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